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The decay rate and the classical radiation power of an 
excited molecule (atom) located in the center of a dispersive 
and absorbing dielectric sphere taken as a simple model of a 
cavity are calculated adopting the Onsager model for the local 
field. The local-field correction factor to the external (radia- 
tion and absorption) power loss of the molecule is found to be 
|3e(w)/[3e(w) + 1]| 2 , with e(u>) being the dielectric function 
of the sphere. However, local-field corrections to the total de- 
cay rate (power loss) of the molecule are found to be much 
more complex, including those to the decay rate in the in- 
finite cavity medium, as derived very recently by Scheel et 
al. [Rev. A 60, 4094 (1999)], and similiar corrections to the 
cavity-induced decay rate. The results obtained can be cast 
into model-independent forms. This suggests the general re- 
sults for the local-field corrections to the decay rate and to 
the external power loss of a molecule in an absorbing cavity 
valid for molecule positions away from the cavity walls. 

PACS numbers: 42.50.Lc, 42.60.Da, 33.50.-j 



I. INTRODUCTION 

Within macroscopic electrodynamics, the decay rate 
F of an excited molecule (atom) at a position ro in an 
absorbing cavity is given by 



2w 



r = "^2" p /* ' ImG ( r 0' r o; w ) • p/*' 



(i) 



with G(r,ro;w) being the (classical) dyadic Green func- 
tion for the system, pfi the relevant dipole transition ma- 
trix element, and u the transition frequency. This result 
is most simply obtained using the classical theory of the 
molecular (radiative) decay |]| in conjuction with the cor- 
respondence principle and extending it straightforwardly 
to absorbing systems ||. Within QED, however, it is de- 
rived by employing the quantized form of the macroscopic 
field in absorbing systems obtained rather recently ||]4|] 
and using the Fermi golden rule plo or solving Heisen- 
berg's equations of motion for the molecule and the field 
in the Markov approximation fJ\M- 

Splitting the Green function into the translationally 
invariant part G°(r — ro;w) and the scattering part 

G sc (r, r ;w), the decay rate can be generally expressed 
as 



r° + r s 



(2) 



where T° is the decay rate as would be in the infi- 
nite cavity medium and T sc is the corresponding cavity- 
induced decay rate. Upon an appropriate regularization 
of G°(r - r ; w), the rate T° is given by |.§ 



r 



free 



3 e"H / c 
2 |e(w)| 2 \ojR„ 



T)(w) 



where 



e [oj 



[f](oj) + in{u>)Y 



is the dielectric function of the cavity medium, R ri 
effective molecule-medium distance, and 



-l free — 



3hc 3 



(3) 



(4) 



(5) 



is the free-space spontaneous emission (SE) rate. The 
first contribution to T° in Eq. (||) is identified as the 
nonradiative decay rate r° r due to the near-field medi- 
ated transfer of the molecular energy to the surrounding 
medium, whereas the second one is the familiar decay 
rate F° ad due to the radiation losses of the molecule in 
the medium, i.e., the SE rate Tg E fH|l(j]. 

For an optically dense cavity medium, the above re- 
sult has to be improved by accounting for the difference 
between the macroscopic field used in its derivation and 
the actual (local) field with which the molecule interacts. 
Restricting ourselves to low-density cavity media, in our 
previous consideration of the molecular decay B as well 
as of the spontaneous emission spectrum [pTj in an ab- 
sorbing planar cavity, we ignored this difference, implic- 
itly assuming, however, that the decay rate corrected for 
the effect of the local field in an isotropic cavity was of 
the form 



f loc = £[r° + r sc ], 



(6) 



where C is an appropriate generalization of the local-field 
correction factor in the lossless case. Thus, as suggested 
by Barnctt et al. for the decay rate T® oc in the infinite 
cavity |5|||, one would have 



^Lor(e) 



e{uj) + 2 



(7) 



in the Lorentz (virtul-cavity) model |h^,y,y,hJJlJ] and 



£ons(e) 



3e(w) 

2e(uS) + 1 



(8) 



in the Onsager (real-cavity) model |]J,|lJ,|19|,gO[ for the 
local field. In either case, the effect of the local field 
would cancel when dealing with the normalized rate T = 

f loc /f? oc = r/r°. 

The conjecture expressed by Eqs. (||)-(||) nas recently 
been shown to be incorrect for both models for the local 
field, which we indicate in Eq . (o ) by putting the tilde 
sign on T. First, Scheel et al. pip demonstrated that, in 
the Lorentz model, a proper inclusion of the (quantum) 
noise polarization in the local field led to a more com- 
plex w-dependence of the decay rate rf oc in an absorbing 
medium than that given by the simple product £loiT 
[see Eqs. (49) and (50) of Ref. |H|]- Their result is in 
full agreement with T® oc in absorbing dielectrics derived 
by Fleischhauer using a microscopic approach J22j. In 
a subsequent work Scheel et al. 17] showed that the de- 
cay rate r° oc in the Onsager model for the local field was 
also much more complex in the absorption case than that 
given by the product £on s r° [see Eq. J27| ) below]. 

Knowing the decay rate rJ ) oc in an infinite cavity 
(medium), it is natural to seek for the decay rate ri oc 
in a finite cavity, i.e., in view of Eq. (||), for the cavity- 
induced rate T^ c . 



Clearly, the calculation of r?? is most 



straightforwadly performed if one adopts the Onsager 
model for the local field. In this model, one assumes 
the molecule in the center of an empty spherical cavity 
of radius small compared with the transition wavelength. 
In view of Eq. (Q), the problem then reduces to the cal- 
culation of the Green function for the system with the 
Onsager cavity at the source position. Owing to the com- 
bined symmetry of such a composite system, however, 
this represents a difficult task, which for most relevant 
cavity geometries demands approximative or numerical 
methods |2JJ. Therefore, in order to make a (first) step 
towards determination of rf£ , in this work we consider 
a special case where the molecule is located in the cen- 
ter of an absorbing spherical particle taken as a simple 
model of a cavity. High symmetry of this configuration 
enables one to perform a simple and exact calculation of 
r^ c in the Onsager model as well as of the power loss 
W^xt 0I the molecule outside the sphere. In turn, these 
results (may) provide a lead to the corresponding solu- 
tions for a more general configuration. We note that this 
special configuration has recently been considered by Cao 
et al. |Q m their calculation of the local-field corrections 
to the SE rate in a purely dispersive medium. We par- 
tially adopt their approach generalizing it to an absorbing 
multilayered spherical system. It should also be noted 
that owing to the existence of high-Q resonances and, 
accordingly, great ability of enhancing optical processes 
p5[ , dielectric microspheres are very attractive objects 
for cavity QED studies. Thus, modification of the decay 
rate and the radiation intensity of an excited molecule 
(atom) in, or near, a (lossless) microsphere has been the- 
oretically considered in both the weak [p|p7 28 29|| and 
the strong |3(J,|8| molecule-field coupling limit, and exper- 
imental observations of modified fluorescence intensity 
have also been reported 31 3^]. 



Owing to its intuitive clarity and easy visualization, 
we prefer a classical discourse in this work. Therefore, in 
Sec. II we recall the classical approach to the problem 
of the decaying molecule and rederive Eq. (J3j) raising, 
through a plausible argument, the question on the com- 
pleteness of this result for absorbing media. In Sec. Ill 
we first obtain formal results for the decay rate and the 
power loss in the external layer of an excited molecule in 
the center of an absorbing multilayered spherical system. 
We then apply these general results to a medium with 
the Onsager cavity and provide an alternative derivation 
of r[ ) oc to that given in Ref. . Subsequently, by consid- 
ering an absorbing dielectric sphere with and without the 
Onsager cavity, we determine T*£ c and W^ in terms of 
the corresponding quantities for the bare sphere. In Sec 
IV we briefly illustrate the effect of the local field in this 
model on T^ c and the total rate ri oc . The main points 
of this work are summarized in Sec. V and the necessary 
mathematical background is given in Appendices A and 
B. 



II. PRELIMINARIES 

In the classical approach, an excited molecule is rep- 
resented by a point dipole pexp(— iujt) at the molecular 
position ro oscillating with the frequency of the transition 
u). The molecular decay rate T is then related through 
r = W/Tiui to the power 



W = -Imp* -E(ro,ro;w) 



(9) 



lost by the dipole in supporting its own field. Equation 
(0) is then obtained introducing the Green function of 
the system through 



E(r,r ;u) 



-G(r,r ;w) • p, 



(10) 



noting that G(r , r ; oj) is the diagonal dyadic and using 
the correspondence principle to let p — > 2p/i |33| . 

The dipole field (Green function) in a cavity can always 
be written as 

E(r, r ; w) = E°(R; u) + E sc (r, r ; w), (11) 

where M 



E°(R;u;) 



3RR 



-(1-ikR) 



i? :? 



47r Yr/T^ .2 1 RR 

— I6(R) + k 2 — - — 



•pe 



'k-n 



(12) 



with I being the unit dyadic, R = r ro, R = R/R, 

and 



k(u>) = V^H- = k'(Lu) + ik"{uj), 



(13) 



is the field of the dipole as would be in the infinite cav- 
ity (medium) and E sc (r, ro;w) is the component of the 
dipole field scattered from the cavity walls. Owing to the 
singular longitudinal (near-field) component of E°(R;o;) 



E°(R;c) 



3RR- I 

R? 



4tt 



i<j(r) 



(14) 



it is the calculation of T° that represents a difficult step 
in determining the decay rate in absorbing systems. 

To remove the singularity from En(R;w) in the spirit 
of the macroscopic-field approach, one may average this 
component of the dipole field over an appropriately 
chosen spherical volume V m = (47r/3)i?^ around the 
molecule or regularize it by letting <5(R) — > 1/V m 
[Spil. In either case, one finds that 



E°(R; W ) 



1 4tt 



r-*o e 3V m 

Since for the dipole transverse field 

E^_(R;w) = E°(R;w) - E||(R;w) 

we find 

u? 2 
E° ± (R;u J )=ik—-p 



(15) 



(16) 



as R — > 0, this leads to the total classical dipole power 

loss 



W° 



^ 4 |P| 2 
3c 3 



3 s"{uj) ( c 



2 \E{LU)\ 2 Vw-Rr, 



?7(a 



(17) 



and, accordingly, to the decay rate T° given by Eq. (||). 
That the above result oversimplifies the frequency de- 
pendence of the molecular decay rate in absorbing media 
becomes clear if one tries to calculate T° using the general 
relation H] 

ImGi,(r,r ;w) = 

If uj 2 
— I d i s—e // (s,uj)Gii(r,s;uj)G* jl (r ,s;uj) 

and applying it to a homogeneous [e(r, ui) — s(lu)] 
medium to determine ImG?- (i"o , ro ; u>). In view of Eqs. 
(H) and (M) , in a less abstract language this is equivalent 
to using the Poynting's theorem to obtain W°. Thus, by 
calculating the dipole energy flow W^ through a spherical 
surface around the dipole and the energy W® absorbed 
per second in the enclosed volume, we find (see Appendix 
A) 



W u 



^ 4 |P| 2 
3c 3 



w \ { i-ikR c y^\ (—r^e 



-2k"R c 



(18) 



where R c is (formally) the lower limit of the radial inte- 
gration in W®- Owing to their characteristic dependence 
on the dielectric function of the medium, we refer to two 
terms in this equation as the absorption (~ s") and the 
radiation (~ 77) contribution to IT" , respectively. Since 
this notation may associate to W® and W2 as the respec- 
tive origins of these contributions, we stress that both W® 
and W9 are needed to obtain each of them, as is clear 
from the derivation in Appendix A. 

To obtain the dipole power loss from Eq. (fL8|), the 
R c — > limit should eventually be taken. However, con- 
sider R c as a small (R c <C A) but finite cutoff for the 
moment. In this case, expanding W° in powers of ujR c /c, 
we find [to 0(^)] 



W° 



^ 4 |p| 2 



3c 3 

2 



kHI 2 



' : ) 3 + e'M ' 



ljR. 



V {u;)e" {w) + n{uj)e' {uj)} 



ujR, 
+ »?(w) 



(19) 



In addition to a near-field term, which very much resem- 
bles the corresponding term in Eq. (O) , now we have two 
new absorption terms. The most striking is the appear- 
ance of another i? c -free term which, therefore, persists 
even in the R c — > limit. This implies that a corre- 
sponding term must appear in any calculation of W° , 
which indicates that Eq. ( |17| ) is, in this respect, incom- 
plete. 

Taking R c in Eq. ([L8|) as a cutoff is equivalent to set- 
ting E°(R;w) = for R < R c . Since this field does not 
obey Maxwell's equations, it is clearly incorrect to re- 
gard W° as the dipole power loss, which we have empha- 
sized by using the tilde. On the other hand, extending 
the macroscopic field down to intcrmolccular distances 
(R c — ► 0) is not justified as the actual field acting on 
the molecule may largely differ from it. These ambigui- 
ties concerning W° are naturally resolved within an ex- 
act macroscopic-field approach in the following section, 
where we adopt the Onsager (real cavity) model for the 
local field and therefore assume the molecule in the cen- 
ter of an empty spherical cavity with the radius R c small 
compared with the transition wavelength A. Since, in this 
case, the longitudinal component EJj(R;a;) of the dipole 
field does not contribute to the molecular power loss, no 
singularity appears in the theory. 



III. LOCAL-FIELD CORRECTIONS 

Consider an excited molecule (dipole) in the center of 
an A-layered spherical system, as depicted in Fig. 1. In 
this case, the power loss W sc of the molecule may be 
written as 



W s 



-Im p*-Ef(r;a;)|, 



(20) 



where Ef c (r;w) = E| c (r, ro; o;)| ro _,o is the scattered part 
of the dipole field in the central region. The calculation 



of the dipole field in this configuration is outlined in Ap- 
pendixB, assuming, for simplicity, that p = pz. Using 



Eq. (|B6|), we find from Eq. 



that (feo = w/c) 



W s 



WfrooReVeTC'f, W ircc = 



ck$\p\ 



(21) 



where C-[ is the corresponding reflection coefficient. In 
the case of the empty central region [W° = Wf roc ], we 
therefore have for the normalized total decay rate T = 
W/Whee of the molecule 



by keep ing only the radiation field (r 
to Eq. dA2|) 



1/r 2 ) contribution 



dfl 



Vn- 



c kp\p\' 

8tt 






- " : A 



' sin 2 0, 



(24) 



t he radiation power VT^f d is given by the last term in Eq. 
(|A3| ) and the total power loss Wn by Eq. jig). Regarding 
the central sphere as the Onsager cavity, we see that the 
local-field correction factor to Wn (and W^ d ) is given by 
\C^/sn\ 2 for £i = l and in the limit fcoai = fco-Rc *C 1. 



f = 1 + Re [Cf ] 



(22) 



Identifying the central sphere with the Onsager cavity, 
this general result provides a direct way for inclusion of 
the local-field corrections to the decay rate in spherical 
multilayered systems. 




FIG. 1. System considered in this paper. All layers are 
assumed absorbing and desribed by the complex dielectric 
function £j(w). 

Of obvious interest is also the radiation power W^f" d of 
the molecule or, generally, the total power loss Wn i n the 
outer reg ion of the multilayer. By comparing Eqs. (Ala) 
and (|B5| ), we see that the dipole field Ejv(r;a;) in this 
region is the same as the field produced in the infinite 
medium (N) of the dipole (C# = C$ + ): 



Pjv 



£N 



(23) 



Accordingly, provided that we let p — ► pn, £ — > £jvj 
and R c — ► rjv— 1, we can adopt all results concerning 
the dipole power loss in an infinite medium derived in 
Appendix A. For example, with these replacements, the 
angular distribution of radiation dW™ /dfl is obtained 



A. Infinite medium 

To calculate local-field corrections to the decay rate in 
an infinite medium with the dielectric function s{u), we 
consider a two-layered system consisting of this medium 
with the Onsager cavity cut around the origin. Then, 
as was done by Scheel et al. [0, If oc is straightforwardly 
obtained using Eq. (E2h and expanding the reflection 
coefficient C 2 (l,£;i? c ) given by Eq. (B7) in powers of 
P\ = k R c < 1, It is very instructive to rederive this 
result by applying the Poynting's theorem to a spherical 
surface around the molecule enclosing the Onsager cavity 
and thus determine the power loss W® oc of the molecule. 
We note that we have performed such a calculation in 
Appendix A to obtain Eq. (18). Thus, W® oc is given by 
this equation provided that we let p — ■*■ p c ff, where [cf. 
Eq. ©] 

Poff = -C|(l,£;i? c )p, 

£ 

and regard R c as the Onsager cavity radius. 

Expanding the coefficient C|(l,e;i? c ) given by Eq. 



(B8) for small p\ — koR c , we find 



Pcff 



3e 



2£+l 



10£ 2 - 9e - 1 



10(2e+l) 



(k Rc 



■f-^^ 3 



■O[(k R c 



(25) 



Owing to the (k R c ) 3 factor in the first term on the rhs 
of Eq. (18), in this term we must use 



(1 - iV^k R c )e l ^ koR c PcS 



3t 



2£+l 



1 14=- 4- 1 

1 + 5 Re ^rr^) 



-21m- 



3 

£^ 



2£+l 



(k R c ) 3 + O[(k R c ) 4 



(26) 



Since the second term on the rhs of Eq. (y_8|) is a well- 
behaved function of k$R Cl in this term it is sufficient to 
let Pcff = 3£/(2e + l)p. In this way, for the normalized 
decay rate fP oc = W?jW hcc we obtain [to 0{*f-)] 



1 loc — 



3e(w) 



2e(u>) + 1 



r?(w) 



kMI 2 



(— ) s 



28|e(w)| 2 + 16e'(tj) + l c 



Hr) 



5|2e(u;) + l| 2 v wi? t 

2 K ( U )|£(w)| 2 + KH£ , (u)+i)He"H 



|2e(t, 



II 



(27) 



The above result coincides with that of Scheel et al. 
pl[ . This time, however, the origin of various contribu- 
tions to Ij c along with their separate local-field correc- 
tions can be clearly identified. As seen from comparison 
with Eq. ( |l7| ) , while the near-field and the radiation-field 
terms get multiplied by £ons as expected || , it is the ap- 
pearance of the additional absorption terms already sig- 
naled in Eq. (n9h that represents essentially new correc- 
tions to r° in absorbing media. Of these two terms, the 
most interesting is the (usually) negative i? c -free contri- 
bution. This term effectively adds to the radiation-field 
contribution (77) to the decay rate in absorbing media and 
tends to diminish the overall rate. 



B. Cavity 

Having determined T® oc , the next step is to consider 
the decay rate ri oc in the general case when the molecule 
is embedded in an inhomogeneous system, i.e., in a cav- 
ity. In order to determine ri oc , we consider the decaying 
molecule in the center of a dielectric sphere of radius R 
and the dielectric function e(u) immersed in an exter- 
nal medium with the dielectric function £ ext (o;). In this 
case, the relevant reflection coefficien t to be inserted in 
11 is Cf{l,e,e cxt ;R c ,R) [Eq. (plCf)] correspond- 



Eq. 

ing to the sphere with the Onsager cavity. Expanding 

this coefficient in powers of pn = koR c , we find that 

C x (1, e, £ext! Rci R) = 
. 9e ,, R s-3 . 9e(8e+l) 1 



2e + l 

Sei ft - 2 
(2e + l) 2 X + 2 



'5(2e+l) 2 
l + O(k R c ), 



(28) 



with /3j given by Eq. (B14). One may recognize that 



2/3i 



01+^2 



-C 1 (e, £ C xt; R), 



where C 2 (e, e cxt ; R) [Eq. (B7)] is the reflection coefficient 
of the system without the Onsager cavity. One may also 
see that 



Re 



9£3 



(2e+l) 5 



35 



2e + l 
18e" 

|2e + l| 



V 

■[(2\e\ 2 + e')K + e" V } (29) 



is equal to the i? c -free contribution to f f oc [Eq. (p7|)] and 
that analogous results hold for the real parts of the first 
two terms in Eq. (|2q). Therefore, from Eq. (E2L) we find 



rioc — i 



loc 



loci 



where 



TL C c = Re 



9ez 



(2s +1) 



;Ci (e, £ C xt; R) 



(30) 



(31) 



is the normalized cavity-induced decay rate with the 
local-field corrections. 

The above result for F^ can be transformed into 



for r° oc . Letting Jl 



a form similiar to Eq. (t 

\JeC\{e, £ cxt ; R) in Eq. (129), we see that Tf^ c is obtained 

from the i? c -free contribution to Ij oc upon replacements 



and k - 

F s 



2A SC , where 
= ReyiiCi (e, £ xt; R) 



(32) 



is the normalized cavity-induced decay rate [Eq. (|2 
and 



A s 



1 



Im7eTC 2 (£,£ext;#) 



(33) 



is the normalized classical cavity-induced level shift |3J] 
of the molecule in the bare sphere. Accordingly, from Eq. 
(E7J) we finally have 



pSC 

1 loc ~ 



3e(w) 



2 i^x 



2e(w) + 1 
2 [2|e(w)| 2 + e'{u)] A sc + e"{u)t s 
|2e(w) + l| 2 



(34) 



We note that this is not an unexpected result once we 
have learned the correct form of the decay rate F° oc in 
the infinite medium. In its derivation it is implicitly as- 
sumed that the radius of the sphere, i.e., the molecule- 
mirror distance, is much larger than the transition wave- 
length (koR 3> 1). Under these circumstances, the 
molecule-mirror interaction goes through the radiation- 
field component of the scattered field and therefore only 
this field component determines the cavity-induced rate. 
Thus, the rate £onsr sc , as would be obtained by letting 
p — > 3e/(2e + l)p in Eq. (pi|), is corrected for an absorp- 
tion contribution in the same way as is the radiation-field 
contribution to the rate F° oc in the infinite cavity. 



C. External region 

To find the local-field corrections to the molecu- 
lar losses in the external region, we consider the field 
Ej,°j(r;w) outside the sphere with the Onsager cavity. 
As already noted, this field is equal to the field in the 



infinite external medium [Eq. ( |Ala ), with k = fc cx t] of 
the dipole [cf. Eq. $2$ " 



loc 
Pcxt 



C 3 (1, s, e cxt ; -R c , -R)p, 

Eext 



(35) 



with the coefficient C|(l, e, e cx t;i? c , i?) given by Eq. 
(pl^). For small p n = k R c , we find that [to O[(k R c ) 2 ]} 



C 3 (1, e, £ ex t; i?ci -R) 



3e 2 



^/efcoi? 2e + 1 /3i + /3 2 ' 



with /3y given by Eq. (B14). Now 



«£oxt 



VefcoR /3l + 02 



-C 2 (e,e cx t;i?), 



(36) 



(37) 



where Cf (e, e cx t; i?) is the external-field coefficient of the 
system without the Onsager cavity [cf. Eq. M, To 
the leading term in kgR c , we therefore have 



loc 
Pcxt 



3<r 



2e + 1 e ext 



C 2 (e,e cxt ;i?)p 



3t 



2e + l 



Pcxt, (38) 



where p GX t is the corresponding effective dipole moment 
for the bare sphere [cf. Eq. (p3|)]. Accordingly, for the 
external field we have 



E^(r;o;) 



3e(w) 
2e(w) + 1 



E ext (r;w), 



(39) 



with E cx t(r;w) being the external field for the bare 
sphere. Consequently, this implies £ons as the proper 
local-field correction factor to the external power loss 



ir, 



loc 



3e(w) 



2e(w) + 1 



W c: 



(40) 



Note that here W ex t is generally given by Eq. (|l8|), with 
fe = fccxt and R c = R. 

We end this section with a remark on a (plausible) 
generalization of the obtained results. Using Eqs. Jig ) 
and (]B6j), the normalized cavity- induced decay rate [Eq. 

can be written in 



(p2|)] and frequency shift [Eq. 
;n function foi 

-p-ImG sc (r ,r ;cj) • p, 



terms of the Green function for the bare sphere as 
3 



pSC 



A s ' 



2fcn 



_3_ 

ik ( 



p-ReG sc (r ,r ;w) • p, 



(41) 



where p gives the direction of the transition. With this 
inserted in Eqs. (fMj), any reference to the specific system 
considered in the derivation of this equation is formally 
lost. Therefore, at first sight, there is no reason why Eqs. 
( p4| ) and (|4l]) should not be taken as the general result 
for the molecular decay rate in an absorbing cavity valid 
for all positions of the molecule away from the cavity 
walls and for all transition dipole orientations. Similarly, 
in view of Eq. (RQ), it appears that Eq. (R9) indicates 



a rather general relationship between the Green function 
elements for a system with and without the Onsager cav- 
ity: 



G loc (r,r ;w) 



3eH 
2e{uj) + 1 



G sc (r,ro;c), 



(42) 



where r is in the external layer and ro in the cavity. If this 
relation holds, £ons(e) as the local-field correction factor 
for the molecular external power loss in the general case 
would then be its immediate consequence. Of course, 
as already stressed, these conjectures cannot be proved 
without the calculation of the exact Green function for 
the system including the Onsager cavity at the source 
position. 



IV. DISCUSSION 

To illustrate the effect of the local field on the decay 
rate in absorbing cavities, we exploit the above simple 
model and consider the decaying molecule in the center of 
an absorbing dielectric sphere surrounded by air (e cx t = 
1). The dielectric function of the sphere is modeled as 



e(u) = e b + 



n 2 



ILOJ 



(43) 



where £(, is the background (high-frequency) dielectric 
constant, ujq and 7 are, respectively, the center fre- 
quency and the width of the absorption resonance and its 
strength is controlled by il 2 . In this analysis, a relatively 
large background dielectric constant Eb — 5 is chosen to 
strengthen the cavity effect of the sphere. 
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FIG. 2. Cavity-induced decay rate with (upper curve) and 
without (lower curve) local-field corrections. The radius of the 
sphere is R = Xo/ir and the material parameters are ej = 5, 
Q, = 0.5u;o and 7 = O.lwo- 

The effects of the local field on the cavity-induced rate 
are illustrated in Fig. 2, where the rates ff^ c [Eq.(pl|)1 



and r sc [Eq.(|32])] are compared for the system with (up- 
per line) and without (lower line) the Onsager cavity, re- 
spectively. The radius (Rlo /c = 2) of the sphere (nearly) 
corresponds to the first peak of T^ considered as a func- 
tion of ujR/c in the lossless (f2 = 0) case and for ui = uj . 
For these parameters, the sphere therefore acts as an en- 
hancement cavity with respect to the SE rate at lo = ujq 
and, since the molecule decays only radiatively, the decay 
rate is equal to the SE rate. Away from the resonance, 
the system is lossless and, owing to the cavity effect of 
the sphere, the two (SE) rates exhibit the familiar oscilla- 
tions with co with the amplitudes of oscillations scaled by 
the usual local- field correction factor Co ns (eb) — 1-85 [cf. 
Eq. (p4|)l . A somewhat different effect of the local field is 
observed in the region of the resonance w ~ u>q. For the 
medium parameters chosen, the second (absorption) term 
on the rhs of Eq. (p4f) is small compared with the first 
one. Moreover, according to our calculations, r^ c cannot 
be distinguished on this scale from £on s r sc over a wide 
range of the parameters fi and 7. The cavity-induced 
rate in this frequency region is therefore (again) predom- 
inantly determined by the product of Cons and T sc . Each 
of these quantities exhibits a characteristic asymmetric 
dispersion around loq. In addition, the sphere is (nearly) 
in resonance with the medium. As a consequence, in- 
stead of a peak at lj = loq appearing in the nonabsorbing 
(£1 = 0) case, owing to the resonant absorption, the rates 
rfo C and r sc exhibit asymmetric double-peak structures, 
which very much resemble the " cavity-polariton" part of 
the SE spectrum in the bad-cavity case [O. 
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FIG. 3. Total decay rate with (upper solid line) and with- 
out (lower solid line) local-field corrections for the system de- 
scribed in Fig. 2. The dotted line represents the macroscopic 
decay rate with the usual local-field correction factor. The 
Onsager cavity radius is R c = 0.1A and R m — R c . 



Figure 3 represents the total decay rates ri oc [Eq. (BC ) 
in conjuction with Eqs. (E7|) and (plj) ] and T [Eq. (|2) 
in conjuction with Eqs. (ft) and (J32J)] with and with- 
out the local-field corrections, respectively, along with 

the decay rate fi oc = £on S f [Eqs. (0) and (0)] with 



the local-field correction factor assumed in our previous 
work |2|,[ll|. In plotting these curves, we have taken that 
the molecule-medium distance R m figuring in T° [Eq. 
(0)] is equal to the Onsager cavity radius, R m = R c . 
Away from the resonance, ri oc and T merely reproduce 
the corresponding curves of Fig. 2, with values en- 
hanced by the off-resonance infinite-cavity decay rates 
f° oc = C 0ns (e b )V^ - 4- 1 and f° = ^ ~ 2.2, respec- 
tively. As before, the situation around coq is different. 
For this R c and the medium parameters chosen, the ab- 
sorption contributions to f i° oc [Eq. @] and f f£ c [Eq. 
(pi|)] are (still) small compared with the corresponding 
radiation contributions. Accordingly, 



floc«£o„ S (^ + f SC ) 



(44) 



holds. For R m ss R c , the same conclusion applies to the 
decay rate T as well, so that 



?/ + r s 



(45) 



The different shapes of the curves in Figs. 3 when com- 
pared with those in Fig. 2 in this region are therefore 
caused by the (anomalous) dispersion of the superim- 
posed Consf] and 77, respectively. This also explains why 
Tioc (upper solid line) practically cannot be distinguished 

on this scale from the decay rate Fi oc = £onsF (dotted 
line), although the corresponding absorption contribu- 
tions to these rates are different [cf. Eqs. (E7J) and (pi)]. 
We note that Eqs. ([hi]) and ( f45| ) hold for a range of 
the Onsager cavity radia and/or medium parameters. In 
these cases, our previous conjecture concerning the local- 
field corrections to the decay rate, as expressed by Eq. 
(0), remains (approximately) valid. 
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FIG. 4. Same as in Fig. 
radius R c — 0.03A. 
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3. but for the Onsager cavity 



With decreasing the Onsager cavity radius, the absorp- 
tion contribution to the decay rate becomes gradually the 
dominant contribution with the leading near-field term 
equal to the nonradiative rate of the molecule. In Fig. 4 
we represent the situation where R c is still large enough, 



so that nonradiative and radiative contributions to ri oc 
are of the same order of magnitude. In this regime, the 
decay rate may approximately be written as 



loc 



-*Ons 



TT?( k oRc)~ 



V 



(46) 



The difference between ri oc (upper solid line) and f 
(lower solid line) this time arises not only because of 
the overall extra factor £ons, but also because of the 3/2 
times larger nonradiative contribution to T [cf. Eqs. ( p7| ) 
and (§)]. This 3/2-factor in the corresponding nonradia- 
tive rates also leads to a significantly larger decay rate 

Tioc (dotted line) than the true rate ri oc (upper solid 
line). We note, however, that the comparison between 

Tioc, T, and ri oc is given here only for illustrative pur- 
poses as their difference is a consequence of our arbitrary 
input: R m = R c . Indeed, in addition to demanding, 
on physical grounds, the same order of magnitude for 
R m and R c , any other relation between these parameters 
may also be assumed. 

Of course, owing to the strong increase of the in- 
dependent terms, for even smaller R c the decay rate 
r° oc largely exceeds the cavity- induced rate T^, so that 
Tioc ~ Tj'oc holds. For a detailed discussion of the decay 
rate in this regime, we therefore refer the reader to Ref. 

0- 



V. SUMMARY 

In this work we have calculated the decay rate and 
the classical external power loss of an excited molecule 
located in the center of an absorbing and dispersive di- 
electric sphere by adopting the Onsager (real cavity) 
model for the local field. We have found that the 
external fields of a dipole classically representing the 
molecule calculated with and without the Onsager cav- 
ity scale with 3e(u!)/[3e(ui) + 1]. This immediately gives 
\3e(co)/[3e(ui) + 1]| 2 as the correct local-field correction 
factor for the (radiation and absorption) power loss of 
the molecule outside the sphere in the absorbing case. 
Whereas this result could have been guessed on the basis 
of a straightforward analytical continuation of the corre- 
sponding result in the lossless case, the local-field correc- 
tions to the total decay rate (power loss) of the molecule 
are found to be much more complex. The total decay 
rate is found to consists of the decay rate for the infi- 
nite sphere, as very recently obtained for an absorbing 
medium in Ref. 0, and of the cavity- induced rate, for 
which we have obtained a similar expression. When ex- 
pressed in terms of the Green function for the sphere, 
these results become formally system-independent. This 
suggests the general results for the local-field corrections 
to the decay rate and to the external power loss of a 
molecule in an absorbing cavity and located away from 
the cavity walls. 



APPENDIX A: DIPOLE POWER LOSS 

To calculate W° using the Poynting's theorem, it is 
convenient to rewrite the dipole electromagnetic field in 
spherical coordinates. With I = ff + 00 + <fi(j> in Eq. 
(O), letting the origin at the dipole site and assuming, 
for simplicity, that p = pz, we have for r ^ (fcp = w/c) 



E°(r;w) = ikk\p- 



' 2h ( i\kr) 
kr 



cosi9f 



^ 1)(fcr) -^(fcr) 



kr 



sm-dO 



B°(r;w) = k 2 k Q ph^\kr) sin #0, 



(Ala) 



(Alb) 



where 



h ( i\z) = -S—, h { l\z) = - e —{l+ % -), 

are the spherical Hankel functions of the first kind. The 
radial component of the Poynting's vector 

P(r;w) = ^-ReE°(r;a/) x [B°(r;u,)]* 



is then easily found to be 



t>, ^ c t l fc l 6 IPl 2 
r • P(r; w) = — Im =±- 



h[ 1] (kr) 



+ h§\kr)] \h^\kr) 



kr 



sin 2 •&. 



(A2) 



This determines the angular distribution of the energy 
flow dWf/dU = r 2 i ■ P(r;w) through a spherical surface 
of radius r around the dipole. Upon the integration over 
the angles, we therefore have 

_ Hp| 2 e" \(l-^kry k r . ^IpI 2 - 2fc -y 

W f\ r )— o UI2 Z?, h? ? o e 



(A3) 

The dipole energy absorbed per second within the vol- 
ume V r of the sphere is given by 



W a °(r) = ^^d 3 r|E°(r; 



-)| 2 - 



(A4) 



With 



|E°(r;-)| 2 = ^ 



6| PP 

2 



h ( ^\kr) 



kr 



h { i\kr) (!) , 
K'{kr, 



kr 



cos 2 d 



sin 2 tf 



and performing the angular integration, we have 



W° a (r) 



w|p| 2 s" 



dr 



Re 



k' 2 -k" 2 -k"\k\ 2 r _ 2k „ 



-j\(l-ikr)e 



ikr I 2 



+ \k\ 4 e 



4„-2fe"r 



(A5) 



Partially integrating the first term and noticing that the 
remaining integrals involving the dipole near field cancel, 
we obtain 



W° a (r) 



cj\p\ 2 e" 



[\ — ikr)e 



ikr 1 2 



^fc 3 lP| 2 (e - 2fc "K c _ e - 2fc " r)- 



(A6) 



This, together with Eq. (A2), gives Eq. (jT|) for the total 
dipole power loss W9(r) + W®{r). 



APPENDIX B: DIPOLE FIELD IN 
MULTILAYERED SPHERICAL MEDIA 

Ge neralizing the result for an infinite medium [Eq. 
(Alb)], the magnetic field of a dipole oscillating along 
the z-axis in the center of an TV-layered spherical system 
[Fig. 1] can be written in the form 



Through 



B(r;w) 



E(r;w) 



eifeop/(r)sma?0. 



V-B(r;w), 



(Bl) 



k e(r) 
the electric field is therefore given by 



£i 



E(r;c)=z^— 
[rf(r)Y 



2/(r) 



;tfr 



r 

sintffl 



(B2) 



with the prime denoting the derivative of the function in 
the brackets. The general form of the function f(r) in 
the Ith layer is 

Mr) = h^ihr^n + C#fc^ \for) + Cf_h^\for), (B3) 

where h\ (z) are spherical Hankel functions. The last 
two terms here are the solutions of the homogeneous 
Maxwell's equations |35[ and give the scattered field 
E sc (r; uj) in the system. The regularity of this field at the 
origin demands that C^_ = C^_ = C± /2, while the out- 
going wave condition at infinity demands that C^_ = 0. 
The rest of the coefficients CK are determined from the 
boundary conditions 



f(rf)=f(rt) and 



[rf{r)]> [rf(r)}> 



z( r i ) 



<r+) 



at the layer interfaces. The scattered field in the central 
layer is therefore generally given by 



2„/~<A r 



E s 1 c {r;u>)=ikik^pCi 



2ji(fcir) 
k\r 



C0S"i?f 



for 



(B4) 



where ji(z) is the spherical Bessel function. In the other 
layers it is given by (E; = E* c ) 



E,(r;w)=tAiftgp 4 /— ^Cj" 



i+ 



2h^\for) 



for 



cos^f 



[hrh^ihr)}' . . 
- sin vu 



fo r 



+Cf 



2h?\k ir ) 



for 



COS$f 



S 2 ), 



for 



(B5) 



with C^„ = 0. Specially, since ji{z)/z — > 1/3 and 
[zjx{z)]'/z — » 2/3 for small z, the scattered field in the 
center of the system is given by 



EJ c (r; W )U = ** 1 ^Cf-p. 

For N = 2, we find (p. L — fori) 

C?(ei,e 2 ;ri) = \ [e 2 h^ ( P2 )[ Pl h^ ( Pl )]' 



(B6) 



eihPMfahPfa)]' 



and 



where 



C| + (ei,e 2 ;ri) 



IE 2 



(B7) 



(B8) 



D = e 1J1 ( Pl )[p 2 h{ 1) (p 2 )r - e 2 h^\p 2 )[p in { Pl )]' . (B9) 

For JV = 3, we have (pij = forj) 

Ci(e 1 ,£ 2 ,e 3 ;r 1 ,r 2 ) = 

'fohPfaJ-fahWjpn) (1) 



Ji(Pn) 



«i/3 2 - a 2 /3i 



-^(Pii) 



C% + (e 1 ,e 2 ,e 3 ;r 1 ,r 2 ) 
C 2 J _(ei,e 2 ,£ 3 ;r-i,r 2 ) 



ft 



ai/?2 - a 2 /3i 

-Pi 
ai/3 2 — a 2 Pi 



, (BIO) 

(Blla) 
(Bllb) 



C|+( £ ii £ 2)£3;ri,r 2 ) 



ie 3 2 

P22 ai/?2 - a 2 Pi ' 



(B12) 



where 



ipu 



«j = — — - [eiji(pii)[p2i/ii (P21)]' 



0), 



-"2 



and 



-e 2 h?(p 2l )[p ll ji(pii)}'], 

(3 j =s 3 h^(p m )[p 22 h[ j \p 22 )}' 
-e 2 h{ j) (p 22 ){p 32 h^(p 32 )}'. 



(B13) 
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